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ABSTRACT 






A branch and exclude algorithm for solution of the 

N N 

"knapsack problem", max £ v,x. where £ w,x. < W and 

0 T 11 II T 11 

1=1 1=1 

x i = 0,1 , is presented which requires relatively small 
amounts of computer running time and core storage allocation. 
In addition, a branch and bound scheme is developed. The 
branch and exclude method is then compared to the branch and 
bound method and to a branch and^baqfld method given by 
Kolesar [2]. Computational results are given. 
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1. 



Introduction 



Stated in precise terms , the knapsack problem is one m 

which one wishes to select from among a finite collection of 

indivisible objects a subcollection which maximizes a linear 

function subject to a linear inequality constraint [2] . 

Dantzig [4] gives the following example: a person is planning 

a hike and has decided not to carry more than 70 pounds of 

different items, such as bed roll, geiger counter, cans of 

food, etc. We try to formulate this in mathematical terms. 

t h 

Let w i be the weight of the l object and v^ be its relative 

value determined by the hiker in comparison with the other 

objects he would like to have on his trip. Let x^ = 1 mean 
t h 

that the i item is selected, and x, = 0 mean that it is 

l 

not selected. We express the weight limitation by 



N 

I w.x. <70, with x. = 0 
. L , i i — i 

i=l 

x. so that the total value 
i 



or 1, and we wish to choose the 
N 



I 

i=l 



v, x. 



i i 



is a maximum. 



The knapsack problem has also been used as a model to 
formulate and solve various aspects of problems arising in 
capital investment [5] , network reliability 16], capital 
budgeting [8], and in optimal methods of cutting stock [7]. 
Dantzig also uses this method in [4] to give a delightful 
"proof" that, of all the forms of marriage (monogamy, bigamy, 
polygamy) monogamy is the best of all possible relations. 

The knapsack problem may be formulated as an integer 
programming problem and solved by any one of the following 
techniques: (i) the cutting plane technique, (ii) techniques 
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employing parallel shifts of the objective function hyper- 
plane, (lii) techniques based upon Boolean algebra, and tiv) 
combinatorial methods which enumerate a restricted subset of 
possible integer solutions. For detailed references on these 
methods see Glover [1] . An alternative formulation of the 
problem may be developed through the use of dynamic program- 
ming f9]. The technique to be presented in this paper falls 
into the class of combinatorial methods mentioned in (iv) 
above . 

In combinatorial methods of this type, the computer 
memory and running time requirements are usually quite large 
for problems involving many variables. 

We present here first a branch and bound algorithm that 
is more in the spirit of Land and Doig [3] than is [2], We 
then develop an algorithm which requires little computer 
storage and running. Computational results are then presented 
for our two algorithms and for Kolesar's method. 
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2. Formulation of the Problem 

We formulate the knapsack problem as follows: 

N 

(1) Maximize J v.x. 

i=l 1 1 



N 

Subject to \ w.x. <_ W, 
i=l 

Where x. = 0, 1 for i=l,2, N 

l 

We can assume without any loss in generality that the 
constants and w^ are integers. Fractions can be handled 
by multiplying through by a proper factor. Non-positive con- 
stants may be handled by the following method which is due 
to Glover [1] : 



i) Let x^ = 



1 

0 



if 


v. > 

1 — 


0 


and 


w . < 

i — 


0 


if 


V. < 

1 — 


0 


and 


w . > 

i — 


0 



ii) Next reformulate the problem as: 

N 

Maximize £ v.y. 

i=l 1 1 

N 

Subject to £ w.y. <_ W , 
i=l 11 



Where y^ = 0,1 for i=l,2, N, 

x . if v . ,w . _> 0 

and where y . = 

1 — x . if v . ,w . <0 

l i i — 

The following example should illustrate this technique: 



Maximize 


3 x i - 


2 x 2 - 4x^ + x 


Subject to 


1 

i — 1 

X 

CM 


3x 2 + x 3 " x 4 



11 



Using Glover's procedure this problem is reformulated 

as : 

Maximize 3y^ + 2y 2 - 1 
Subject to 2y^ + 3y 2 £ 4 

and the optimal solution if the integer constraint on is 
relaxed is 10/3 units with X = (1,1/3, 0,1). 
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3. A Branch and Bound Algorithm 

As shown by Dantzig [4] the optimal fractional solution 
to 

N 

(2) Maximize ) v.x. 

i=l 1 1 

N 

Subject to £ w.x. < W 

i=l — 

0 1 x i < i ( i=l , . . . ,N) 

is given by 

x. = 1 if i < r 

l 

x^ = 0 if 1 > r 
W-£ w . 

* x r = 1 < r 1 

w 

r 

where r is the least integer (0 _< r <_ N) for which 

J w. > W . If no r exists we have all x. = 1. If 

“ l — l 

i <_ r 

x r = 0^ we have the optimal solution to (1). 

If x r is fractional the value of the objective function 

is z(l) = £ v. + v x . We consider z(n) , the value of 

i < r 1 r r 

the objective function at node n, as the solution to (2) with 
assigned variables added as constraints. A branch and bound 
scheme is essentially as follows: 

1. Label node 1 with z(l). Go to 2. 

2(a). Find the terminal node with the largest value of 
z (n) . This is the node at which the next branch- 
ing will take place. Any node 1) contains the 
effect of assigning values to variables and solv- 
ing (2) with the assigned values of the variables 
added as constraints. 



13 



(b) . If the solution at node n has all integer var- 
iables, we have achieved an optimal solution to 
(1). Stop. If not, go to 3. 

3(a). Set n = n+1 and some unassigned variable, say x^, 
equal to 0. Solve (2) with all assigned variables 
added as constraints. Label node n with the value 
z (n) . Go to 3(b). 

(b) . Set n = n+1 and x^. = 1. Solve (2) with all 

assigned variables added as constraints. Label 
node n with the value z (n) . Go to 2. 

Kolesar's algorithm consists in taking as x^ , in step 
3(a), the unassigned variable with smallest index (i.e., 
t = i for which v^/w^ is the maximum for unassigned variables 

V- 

We propose instead to follow [3] and to take x fc , in step 
3(a), the variable that is fractional at node n. 

To illustrate our branch and bound algorithm, we solve 
the problem given in [2]: 



Index 


w . 

l 


V. 

l 


1 


30 


60 


2 


50 


60 


3 


40 


40 


4 


10 


10 


5 


40 


20 


6 


30 


10 


7 


10 


3 



and W = 100. Solving (2) , we obtain z(l) = 140 with x^ = 1 , 
x^ = 1, x 3 = 1/2. We label node 1 with z(l) = 140. Since 
x^ is fractional we branch from node 1 and proceed to step 
3(a). We solve (2) with x^ = 0. This produces z(2) = 135 
with x 1 = 1, x 2 = 1, x 3 = 0 , x 4 = 1, x 5 = 1/4. We label node 
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Figure 1 # Complete tree for the example 



15 




2 with z(2) = 135. In step 3(b) we set x^ = 1 . This produces 
z(3) = 136 with = 1 , X 2 = 3/5, x^ = 1. We go to step 2 in 
the algorithm. We see that z(3) is the maximum for all ter- 
minal nodes. Since the solution at node 3 has X 2 fractional 
we branch from node 3 and proceed to 3(a). We solve (2) with 
x^ = 1 and X 2 = 0 added as constraints. This produces 
z(4) = 120, with x^ = 1 , X 2 = 0, x^ = 1 , x^ = 1, x^ = 1/2. 

We proceed to 3(b) and solve (2) with x^ = 1 and x 2 = 1. 

This produces z(5) = 120 with x^ = 1/3, x 2 = 1 , x^ = 1. Re- 
turning to step 2 we see that z(2) is the maximum for all ter- 
minal nodes. Since the solution at node (2) has x^ fractional 
we must branch from node 2. The method continues easily and 
the complete tree is shown in Figure 1. The optimal solution 
is given at node 8 with z(8) = 133, x^ = 1 , x 2 = 1 , = 1 , 

x^ = 1 and all other x^ = 0 . 

The algorithm achieved solution with the generation of 
9 nodes. The same problem was solved in [2] and required 15 
nodes. Further computational experience is presented in 
Section 5. 
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4. A Branch and Exclude Algorithm 

Branch and bound techniques normally require large 
amounts of computer storage because of the necessity to store 
information for all terminal nodes. In this section, we pre- 
sent a branch and exclude algorithm for solving (1) that re- 
quires little computer storage. 

Our branch and exclude algorithm first finds an obvious 
integer solution to the constraints of (1) . This solution 
is a lower bound to the optimal solution. We develop a branch 
of a tree and explore each part of the branch until the lower 
bound is reached or until a new feasible solution is found 
that represents a larger lower bound. We then back- track 
and develop new branches of the tree developing possibly 
larger lower bounds. Further branching is excluded when the 
lower-bound is reached. The algorithm stops when all new 
branches are excluded. The only information that is stored 
is the current lower bound solution and the branch routing. 

At the end, the lower bound solution is optimal. 

Define [x] as the greatest integer less than or equal 

to x. 

Define S (Xq , x^ , . . . , x^ as the current lower bound solu- 

N 

tion, where the x. are all given and x n = Y v.x. . 

' 1 3 0 . L T 1 1 

1=1 

Define X (x^ ,x 2 , . . . ,x N ) to indicate assigned variables. 

A value x^ = 2 (or any number not equal to zero or one) in- 
dicates the variable is unassigned. An assigned variable 
will have the value zero or one. 

th 

Define R(L) as the index of the L assigned variable. 
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The solution to (2) , with L assigned components of X 



added as constraints, is 

x^=l if i < r, i / R(j) 
x i = 0 if i > r, i ^ R(j) 
(3) X R(j) = 0 or 1 (j = 

assignment 



(D — if • • • f L ) 

( j = 1 , . . . ,L) 
depending on the 



x 



r 



(W 



L 

w R(i) x R(i) 



- I W. )/w 

ieM(L) 



Z (L+l) 



L 

V R(i) X R(i) 



+ l V. + 

ieM(L) 1 



v x 
r r 



where the set M is given by 

M(L) = { i | i < r, i j- R ( j ) (j = 1 , . . . ,L) } 
and r is the least integer (0 < r < N) for which 



L 

l 

ieM(L) 



w . + w 



> w - T 

r — . L 



i=l 



W R(i) X R(i) 



If no r exists we have all x^ = 1 for i / R ( j ) (j 



1 , . . . , L ) . 



A lower bound to the solution of (1) is given by 



x (L) 



L 

J 1 V R (l) X R ( i ) 



+ l 

ieM(L) 



v- 

l 



The algorithm follows: 

1. Set L = 1 and all components of X to two. Go to 2. 

2. Solve (2) . If 'the solution is all integer we have the 

optimal solution. Stop. If not the solution is x^ = 1, 

x 2 = l,...,x ^ = 1 and x f is fractional. We calculate 

Xq = J v. and form S (x^ , 1 , 1 , . . . , 1 , 0 , 0 , 0) as a lower 
i<r 1 

bound to the optimal solution, where the zero components 

4“ V» 

in S represent x^ = 0 , i >_ r. Set R(l) = r and the r 
component of X to zero. Go to 3. 
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3(a). Solve (2) with the L assigned components of X added 
as constraints. 

We obtain (3). If (Z (L + 1)] <: x Q to to 4. If 
Z(L + 1) > Xq and we have an integer solution, take 
x 0 = + 1)/ form a new S (Xq ,x^ , . . . ,x^) from (3), and 

go to 4. If [Z (L + 1) ] > Xq and we do not have an integer 
solution, go to 3(b). 

(b) . If x(L) > Xq , take Xq = x(L) and form a new 
S ( x 0 ,X 1 , * * ‘ ,x n) f roin (3) with x r = 0. In any case. Set 
L = L + 1, take R(L) = r, and set the R(L) component of X to 
zero . Go to 3 (a) . 

4(a). If the R(L) component of X is equal to zero change the 
component to one and go to 3(a) . If the R(L) component of 
X is one go to 4 (b) . 

(b) . If L = 1 the optimal solution is S (xq , x^ , . . . , x^) . 

Stop. If L ^ 1 change the R(L) component of X to two, set 
L = L-l , and go to 4(a). 

This completes the algorithm, z (L) represents the value 

th 

of the objective function at the L level of a branch. R(L) 
represents a routing of assignments along the branch. Only 
one branch is studied at a time with preference given to 
assigning the value of zero to the fractional variable. This 
allows lower-bounds to be achieved more rapidly. The only 
permanent storage information required is the current lower 
bound solution S, the assigned variables R(L) at level L of 
the branch, and the assignment vector X. 
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In the problem given above the algorithm as viewed in 
Figure 1, goes as follows: 

1. At node one the optimal solution to (2) is 
fractional. Xq = 120 . 

2. Branch to node two, assigning x^ = 0. 

3. Branch to node six, assigning Xg = 0 . 

4. Branch to node eight, assigning Xg = 0 . This 

produces an all integer result with Xq = 133, 
x x = 1, x 2 = 1, x 4 = 1, x ? = 1 . 

5. Back-track to node six; branch to node 9, 
assigning Xg = 1. 

6. Back-track to node 6, removing the assign- 
ment on x, . 

b 

7. Back-track to node 2; branch to node 7, 
assigning Xg = 1. 

8. Back-track to node 2, removing the assign- 
ment on Xg . 

9. Back-track to node 1; branch to node 3, assign- 
ing Xg = 1. Calculate x(l) = 100; do not 
change the current lower bound. 

10. Branch to node 4, assigning Xg = 0. 

11. Back-track to node 3; branch to node 5, assign- 
ing Xg = 1 . 

12. Back-track to node 3, removing the assignment 
on x 2 . 

13. Back-track to node 1; end. 

14. The optimal solution is given at node 8. 

In this example, the total tree investigated has the same 



number of nodes. 
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5. Computational Results 

We have programmed the two branch and bound schemes and 
the branch and exclude algorithm in Fortran IV and have run 
test problems on the IBM 360/67. A program for the branch 
and exclude method is given in Appendix A. 

The coefficients v and w^ were generated as random 
integers and the three methods were then tried on the same 
problem. If 

BK = Kolesar's branch and bound algorithm, 

BB = the branch and bound algorithm presented here, 

BE = the branch and exclude algorithm presented here, 
we obtained the following results for the number of nodes 
generated and the time to achieve solution: 



TABLE 1. Comparison of Methods 



w 


N 


20 


30 


40 


50 


20 


30 


40 


50 




BE 


28. 8 


28.0 


28.8 


35.7 


61 


87 


120 


164 


25 


BB 


26.4 


20.6 


26.6 


28.6 


201 


282 


395 


494 




BK 


56.2 


49.6 


1 

loo 

• 

IO 


Ol 

• 

Ol 

in 

i 


259 


345 


447 


552 




BE 


34.5 


38.7 


36.1 


54.5 


75 


111 


145 


235 


50 


BB 


38.6 


48.2 


- 42.0 


51.8 


236 


361 


453 


581 




BK 


62.2 


78.4 


78.0 


92.6 


276 


434 


553 


723 




BE 


19. 8 


39.0 


27.8 


18.5 


58 


116 


128 


135 


75 


BB 


31.0 


30.2 


27.2 


25.0 


222 


312 


404 


485 




BK 


62.6 


101.2 


84.4 


61.6 


280 


523 


578 


616 




BE 


7.2 


24.3 


46.8 


29.0 


36 


97 


181 


163 


o 

o 


BB 


16.6 


24.8 


64.0 


50.4 


191 


301 


560 


598 




BK 


56.6 


65.4 


136.6 


134.2 


273 


40 3 ' 


806 


978 



Average number of nodes generated. Average time in Milli- 
seconds. (10 problems run to obtain each average) 
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In these test problems superiority of BB and BE over BK 
is clearly apparent with regard to number of nodes generated. 
As for time considerations, the branch and exclude^ method 
shows to be markedly faster over the other two methods, while 
BB is faster than BK. Various other problems than those 
listed above were run, and in general, the BE method achieved 
solution in about one-half to one-third the time required by 
the BK and BB methods. 

The branch and exclude method was tried on several prob- 
lems with 5000 variables and solution was achieved in all 
cases in approximately 4 minutes. The other two methods 
failed because of storage limitations. In general, the branch 
and exclude algorithm presented here can rapidly solve ex- 
tremely large problems. 
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APPENDIX I 



FORTRAN IV Computer’ Program for the Branc 1 ' and Exclude 

Method 



LEVEL 1* MOO 0 



MAIN 



DATE * 68051 



C 

C 

c 

c 

c 

c 

c 

c 

c 



IMPLICIT REAL! Y-Z,i) , INTEGER I A- X I 

INTEG€ft*2 V.W.OROER.X, SVECT.H.F ,G,H 

COMMON VI 1001 ,W< 1001 .N, XI 100 I , OROER 1 100 ) , WNOT , 

♦YVEC I 100 I.SFLAG.XNOT.SVECTI 100 ) , L ,R 1 100 ) , FI 100 ) . 

*G I 100 I 

DIMENSION HI 100 I 

V » W » R » X ANO L ARE AS STATED IN THE TEXTUAL EXPLANATION OF 
THIS ALGORITHM 

WNOT IS THE WEIGHT CONSTRAINT. W 
N IS THE NUMBER OF VARIABLES 

SVECT IS THE SOLUTION VECTOR. S. EXCEPT THAT X-SUB-ZERO 
IS KEPT AS XNOT 

SFLAG IS THE NUMBER ASSIGNEO TO A PARTICULAR NOOE 



READ 101.WNOT.N 

101 FORMAT (2151 
00 102 1 = 1. N 

reao loi.vm.wm 

102 CONTINUE 

CO 103 1=1. N 
R C I 1=0 

xm=2 
SVECTl I 1=0 

103 CONTINUE 
XNGT=0 
SFLAG=0 
L* 1 



C 

C THE VARIABLES ARE HERE PLACED IN ORDER OF DESCENDING 
C VSUB I/WSUBI 

C 

CALL RANK 
C 

C THE OPTIMAL SOLUTION IS NOW FOUND BY THE SOLVIT ROUTINE 
C 

CALL SOLVIT 
C 

C AT THIS POINT THE VARIABLES ARE REARRANGED INTO THEIR 
C ORIGINAL ORDER. AND THE RESULTS PRINTED 

C 

DO 91 M.N 
F I CRDER ( 1 1 l = SVECT< II 
G I CROER ( I ) ) = VI I ) 

HI CRDER ( I ) l=W< I ) 

91 CONTINUE 
PRINT 70, XNOT 

70 F0RMATUH1, 20X, 'SOLUTION =*, 15,' UNITS' ) 

PRINT 90 

90 F O RMA T ! 2GX f » I • r!9X, »X< I » WHrX . ' VI It • , L6X, ' W< I M ) 

00 92 1=1, N 

PRINT 93, I, Fill, G( II, HI It 
93 F0RMATI3X.AI17X.mi 

92 CONTINUE 
STOP 
END 
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nocooo v\ * non nj oooh noon 



i V-tr* ) » IN TE G ER I A - X 

I, ORDER » , " “ 



INTEGER*2 V ,W, ORDER, X, SVECT,R,F ,G,H 
COMMON V 1 100 I . W ( 100 I , N , X ( 100 ) , 0R0E R < 100 I • WNOT , 

♦YVECI 100 ),SFLAGtXN0T»SVECTI 100 ) ♦ L, R ( 100 1 , F I 100 1 , 

*G( 100 ) 

ORDER IS THE VECTOR WHICH WILL CONTAIN THE ORDER IN WHICH 
THE VARIABLES ARE REARRANGED 

DO l 1*1, N 
GROER <11*1 
Y*WII) 

z*vm 

YVEC( 1 1 *Z/Y 
CONTINUE 

THE RANK GRO E R IS NOW DETERMINED AND STORED IN ORDER 

K=0 
K = K+1 
J*K+ 1 

ZAP*YVEC(K) 

M*K 

DO 3 I* J , N 

TFI-V V E CI I T, LE« ZAP) GO TO 3 
ZAP=YVEC(I ) 

M* I 

CONTINUE 

YVECIN)=YVECIK) 

YVEC ( K l*ZAP 
PI P=0RDER(M) 

CROERIM)=ORDERIK) 

CRD E R IK )= P I P 
IFIJ,NE.N)GO TO 2 

THE VARIABLES ARE NOW REARRANGED 

DO 4 1*1, N 

FU)*vm 

gui=wu> 

r n KiT I |U| I P __ 

Lutr i t nuc 

DO 5 1*1, N 

VI II=F|0R0ERI I) ) 

MCI »=G(0RDER( 1)1 

CONTINUE 

RETURN 

ENO 



SUBROUTINE SOLVIT 

IMPLICIT R E AL T V -Z^H-vENTEGEIHA-X ) 

I NTEGER*2 V, W, ORDER, X,SVECT,R,F ,G,H 
COMMON VI 100),WI 100) ,N,X| 100 ), ORDER ( 100 ) , WNOT , 
♦YVECI 100 ) , SFLAG, XNOT, SVECTI 100) ,L,R( 100) ,F( 100) , 
*G< 1001 

SUM IS THE VALUE OF THE SUM OF WSUBI*XSUBI FOR ALL I 
TOTALIS THE VALUE OF THE SUM OF VSUBI*XSUBI FOR ALL I 

Y ^WttK", S" igi« 18 Vg 8 ' 
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10 

11 

3 



C 

C 

C 

C 



«WI,W’ N 

CONTINUE 
SUM=0 
T0TAL=0 
SFLAG-SFLAG+l 

THE VARIABLES WHICH HAVE BEEN CONSTRAINED TO BE ONE ARE 
HERE CONSIDERED 



00 12 I 



I = 1 » N 
IFIXI n.NE.lJGO TO 

if)*! 



12 



12 

C 

C 

C 

C 



YVECI 

SUM -SU M+W I I | 

TCTAL=TOTAL+Vt 1 1 

IF( SUM-WNOT 112,14,20 

CONTINUE 

IF THE CONSTRAINT EQUATION HAS NOT BEEN VIOLATED* 
UNASSIGNEO VARIABLES ARE NOW CONSIDERED 



C 

C 

C 

C 



00 13 1 = 1, N 

IFI X I D.LE.IIGO TO 

yveci n=i 
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THE INITIAL LOWER BOUNO VECTOR, S, IS BEING FORMED BY THE 
NEXT CARD 

IFI SFLAG.EQ. 1) SVECTI I 1 = 1 
SUM-SUM+WI n 
TGTAL=TOTAL+V(II 
IFI SUM-WNOTI 13,14, 15 
CONTINUE 



13 
C 

C IF 

C 

C 

14 



STATEMENT 14 IS REACHEO, AN INTEGER SOLUTION HAS BEEN 
FOUND 

CAttr~ SYORGI TOTAL) 

IFI SFLAG«EQ« 1 1 RETURN 



C THE VALUE ASSIGNED TO THE FRACTIONAL VARIABLE IS NOW 



C 

C 

16 

20 



17 



CHECKED, ANO BRANCHING INSTRUCTIONS GENERATED 
G=X( R I L ) ) 

IFIQ.EQ.OIGO TO 17 
X ( R ILII »2 
L=L- 1 

IFIL.EQ. 1 (RETURN 
GO TO 16 
XIRILI 1=1 
GO TO 10 



C IF STATEMENT 15 IS REACHEO, A NON-INTEGER SOLUTION HAS 

C Bt€N F OU N D 

C 

15 SUM*SUM-Win 

TOTAL=TOTAL-V( 1 1 
YVECI ll»0 

C 

C THE NEXT CARO CHECKS TO DETERMINE IF A NEW LOWER BOUND 
C CAN BE ENTERED. 

C 1 F I M R R » I « ED, 1 . A NB . TOTAL »GT . X N OT I CALL STO R E ♦ TOT AL T 

C 

C FORMULATION OF THE INITIAL LOWER BOUNO IS COMPLETED BY 
C THE NEXT THREE CAROS 

C 
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nonnnno wt-**-* noooooo^ n> •goonoon 




FLAG.NE. 1)60 TO 79 

:n*o 

TOTAL 



9 



A DETERMINATION IS NOW MADE AS TO WHETHER THE GREATEST 
INTEGER CONTAINED IN THE VALUE OF THE PRESENT 
NON-INTEGER SOLUTION IS EQUAL TO THE VALUE OF THE 
LOWER BOUND SOLUTION 

Z*WNOT-SUM 



v»wm 

vv e cc 



ff-Z/Y 

IFCSFLAG.NE.DGO TO 2 
ZSOL*TOTAL+VC I )*YVECC I ) 

SI NT=TOTAL + C WNOT-SUM ) *VC I }/M(II 
IFCSINT.EQ. XNOT IRE TURN 
GO TO 6 

Z*TOTAL+VCI)*YVECCI I 
SINT*TOTAL+CWNOT-SUM)*VC I) /W (I) 

IFC SIN T.Lfr^XNOT LGO TO 1G - 

IF THIS POINT HAS BEEN REACHED* THE LARGEST INTEGER 
CCNTAINEO IN THE VALUE OF THE CURRENT NON-INTEGER 
SOLUTION IS LARGER THAN THE VALUE OF THE LOWER BOUND 
SOLUTION* AND SO THE FRACTIONAL VARIABLE IS ASSIGNED 
A ZERO VALUE 



OQ -I B- J~1 *N 

IF<YVEC(JI.LT. I. AND. YVECC J ) .GT.O )G0 TO 

8 CONTINUE 

9 L*L+1 
CONTINUE 
R CL )*J 
XC J)*0 
GO TO 10 

END 
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«■ 

INTEGER*2 V.W,OROER, X* SVECT,R*F,G,H 
COMMON VC 100) *WC 100) ,N,XC 100 ) .ORDER C 100 ) , WNOT* 
♦YVECC100) , SFLAG* XNOT * SVECTC 100 )*L*RC100)*FC100), 

*GC 100) 

IF THIS ROUTINE HAS BEEN CALLED. TOTAL IS THE VALUE OF AN 
I NTEGE RESOLUTIO N 

CURRENT LOWER BOUND SOLUTION, AND THE SMALLER OF THE 
TWO SOLUTIONS WILL BE REPLACED BY THE LARGER 

IF C TOTAL #LE^ XNOT) RE TURN 
XNOT-TOTAL 
00 4 1*1, N 
SVECT C I )*YVEC( I ) 

CO N TI N U E 



RETURN 
ENO 
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APPENDIX II 



Logic Diagram for the SOLVIT Subroutine Contained 

in Appendix I 
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